O.P.E. and Power Corrections to the QCD coupling constant. 
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Lattice data seems to show that power corrections should be convoked to describe appropriately the transition 
of the QCD coupling constant running from U.V. to I.R. domains. Lhose power corrections for the Landau- 
gauge MOM coupling constant in a pure Yang-Mills theory (Nf = 0) are analysed in terms of Operator Product 
Expansion (O.P.E.) of two- and three-point Green functions, the gluon condensate < A 2 > emerging from this 
study. Lhe semi-classical picture given by instantons can be also used to look for into the nature of the power 
corrections and gluon condensate. 



1. The non-perturbative QCD coupling 
constant from lattice QCD. 

The non-perturbative calculation of the QCD 
coupling constant is a major open question, as 
it could give some insight into the problem of 
confinement. This calculation, as a very impor- 
tant task, has already been performed by sev- 
eral different methods, The method used 
here is based in the calculation of the coupling 
constant from the three gluon coupling ]^), that 
allows an easier physical interpretation than the 
Schrodinger functional method, for example, and 
if it does not allow to study a quite wide range of 
energies, it will allow us to focus on the regime of 
energies where the transition from perturbative 
to non-perturbative regimes in QCD occurs. 

The calculation of the coupling constant in 
this method comes from the lattice evaluation of 
two- and three-point gluon Green functions, the 
renormalised coupling being defined (according to 
MOM schemes) as the renormalised three gluon 
vertex where external propagators are explicitly 
amputated: 
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9(k 2 ) = 



GW(p\)GW(pi)GV)(pi) 



(1) 



where G^ 2 \p\) and G^\p\,p\,p\) are two- and 
three- point Green Function form factors |i|,Q, 
and Zz(k 2 ) = k 2 G^\k 2 ) is the gluon propagator 
rcnormalization constant. 



2. O.P.E. and (A 2 ) condensate. 

Lattice calculations of the QCD coupling con- 
stant suggest the necessity to add power correc- 
tions to the purely perturbative expressions, to 
correctly describe its running An Operator 
Product Expansion (O.P.E.) analysis of the Green 
functions in Landau gauge^j relates this power cor- 
rections to the existence of a non-perturbative 
(A 2 } condensate ||, through expressions: 



1 In the lattice we will work in the minimum A 2 Landau 
gauge, c9 M A M = (Absolute Landau Gauge), so all gauge 
dependent quantities will be expressed in this particular 
gauge. 
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Go^.P.eXP 2 ) — Gplrt.iP 2 ) + C ^2 ' 

af- p - E (p 2 )^a^(p 2 )+d { ^-^, (2) 

where perturbative expressions are developed at 
three loops, and the functions c and d include 
the Wilson coefficient of the expansion and the 
anomalous dimension of the condensate at leading 
logarithm. 

By performing a combined fit of lattice results 
in a wide region of energies (from 3 to 10 GeV) to 
expressions in (|J), in two different MOM schemes, 
a value of A^jg is extracted, in fairly good agree- 
ment with the one obtained by the ALPHA col- 
laboration Q, by a completely different method. 
A value of the {A 2 ) condensate comes out from 
the analysis. 

The physical meaning of this condensate is still 
an open question, and a lot of work is being de- 
voted to its study during last years, for exam- 
ple, in relation to confinement The aim in 
this work will be to study the possible semiclassi- 
cal contribution to this condensate coming from 
instantons, and whether they might explain the 
presence of power corrections in Green Functions. 



A« = (AJ) Jnat . + QJ((A-) /nat .) , (3) 

where in principle quantum fluctuations could 
deppend on the semi-classical background. 

If we nevertheless assume the hypothesis that 
quantum fluctuations are not sensitive to the 
background ||, the O.P.E. (A 2 )-condensate, cor- 
responding to the non-perturbative part of (A 2 ), 
will be: 

( a2 )o.p.e. « WLst.) « (A 2 Q 2 (0)) . (4) 

The task in this work will be to compute the 
instantonic contribution to this condensate, that 
(up to other non-perturbative contributions, that 
could of course exist) we will finally compare to 
the O.P.E. one, computed in Q. 

In order to compute (A 2 ) in an instantonic 
background, we will consider the gauge field in 
an ensemble of non- interacting SU(3) instantons 
(I) and anti-instantons (I) with radii {pi}, cen- 
tered in {zi}, and with an orientation in color 
space given by {{O aa )i}. If they are not inter- 
acting and randomly oriented, their contribution 
to the (A 2 ) condensate in the absolute Landau 
gauge would give: 



3. The role of instantons. 

Instantons are classical solutions of QCD 
equations of motion that have been longly stud- 
ied as a possible description of QCD vacuum, and 
claimed to explain some non-perturbative proper- 
ties of QCD at low energies, as the axial anomaly, 
rj mass, chiral condensate, etc (See § for a gen- 
eral overview). 

Based in the fact that instantons are saddle 
points of QCD action, we could think of factoriz- 
ing the path integral into an integral over semi- 
classical gauge field configurations (in this case 
instantons) and the integral of quantum fluctua- 
tions around this semi-classical background. This 
means that the gauge field (A^) could be decom- 
posed into a instantonic background ((A®)i ns t.) 
plus quantum fluctuations: 



N f 

(Anst « v J = 12 *V». (5) 

where A"(x) is the standard 't Hooft Polyakov 
instanton gauge field 0, p the average radius, 
and n = N '^ N i t ne density. 

If we accept the phenomenological values as- 
signed to n and p by the Instanton Liquid Model 
H (n ~ 0.5/m~ 4 and p ~ 1/3/m), the instan- 
tonic contribution will be (A 2 ). . ~ 0.5GeV 2 . 

\ / Inst. 

We will perform, however, our own analysis, thus 
testing the latter approach. 

3.1. Cooling. 

In principle, a direct measure of A 2 in the lat- 
tice should be possible, but the presence of the 
UV divergent part is hardly separable from the 
soft, instantonic one. The other possibility is to 
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perform a cooling procedure, that, through a pro- 
gressive elimination of UV fluctuations, will allow 
us to compute the number and size of instantons, 
giving an indirect measure of the {A 2 ) conden- 
sate through (||) . Both methods agree after a high 
number of cooling sweeps, when quantum fluctua- 
tions are almost completely suppressed, but then 
cooling bias is so strong, as we will see, that no 
information can be recovered about the original 
situation. 

We will use the traditional cooling method jlO| , 
even if it introduces a number of known biases, as 
I — I annihilation, and a modification of instanton 
sizes and lattice spacing. The approach proposed 
here is to compute instanton properties for dif- 
ferent number of cooling sweeps, and extrapolate 
back to the thermalised situation |ll| , in order to 
recover their physical meaningP]. 

3.2. Shape Recognition. 

Instantons will be localised in cooled lattices 
via a geometrical method (Described in |l3|] .) 
that accepts a topological charge lump as an in- 
stanton when the ratio of the integral over lattice 
sites with topological charge bigger than a given 
fraction, a, of the topological charge at the max- 
imum and its theoretical counterpart, 



I , \QpM\ > d 4 xQ p (x) 

1 IQ P (Q)I - 

1 - 3a 1 / 2 + 2a 3 / 4 



(0) 



is around 1 for a range of values of a. To verify 
that it is really an instanton, we compute (Q) both 
for the topologial charge density and the action, 
accepting it as an instanton only if it is self-dual. 

Once the lump has been identified as an instan- 
ton, the radius will be computed from the size 
of the cluster where the integral has been devel- 
oped, in agreement with the value of the radius 
deduced from the value of the topological charge 
at the maximum. 

We have tested this method for different lattice 
spacings, from 0.055 to 0.23 fm. As a matter of 
fact, for big lattice spacings, instantons have radii 



2 The usj 
oped in [ 



f improved cooling methods, as the one devel- 
, could improve this approach, as radii evolu- 
tion is minimised, but I — I annihilation is unavoidable, so 
the extrapolation will be anyway necessary. 



of order of the lattice spacing, being in practice 
indistinguishable of UV fluctuations. Therefore, 
only lattices finer than 0.1 fm will be used, oth- 
erwise an important number of instantons with 
small radii will be lost (all for radii below two 
lattice spacings). 

3.3. Results. 

If we compare the measures of the instanton 
density for two different lattice spacings, 0.074 
and 0.101 fm, and different volumes, as a func- 
tion of the number of cooling sweeps (Figure , 
we can conclude that no finite volume effect af- 
fects our results, as the measures for the 24 4 lat- 
tice and 32 4 one are compatible. This is true ex- 
cepting lattice volumes approaching the instanton 
volume. 



• L=24, p=6.0 
o L=32, (i=6.0 
□ L=32, (5=6.2 




Figure 1. Density of instantons measured form dif- 
ferent lattice volumes and lattice spacings (a((3 = 
6.0) = 0.101 fm, a{(3 = 6.2) = 0.074 fm; 



For different values of the lattice spacing, on 
the contrary, bigger densities are found at the 
same number of cooling sweeps for finer lattices. 
In addition to the problems to detect instantons 
with small radii, there is a known problem of cool- 
ing procedures, they are not independent of the 
lattice spacing Jig] , and both effects could play 
an important role to explain the difference found 
between different lattice simulations. 
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To illustrate this, if we fix the value of the den- 
sity (that is, we perform a different number of 
cooling sweeps for each lattice spacing) and com- 
pare the average radii, they agree within the er- 
ror limits. Moreover, comparing the distribution 
of instanton radii, they are in a quite good agree- 
ment (Figure |J), even with very few statistics (10 
configuration each). 
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Figure 2. Comparison on instanton histogram at 
n c = 10 for P = 6.0 and n c = 15 for /3 = 6.2. 



3.4. A naive model of annihilation. 

With the method outlined above, we compute 
the density and size of instantons in a lattice, for 
different numbers of cooling sweeps, n c , obtain- 
ing values with a strong dependence on n c (See 
figures) , that avoids to obtain any physical infor- 
mation at fixed n c . Moreover, this dependence 
on the number of cooling sweeps is different at 
different lattice spacing, so an extrapolation of 
results after cooling is unavoidable to obtain any 
physical information. 

As a first approach to the understanding of this 
evolution, we will make a simple model, where 
instantons annihilate with antiinstantons (Being 
so AN = Ni — Nj a constant) proportionally to 
their packing ratio, and to the number of antiin- 
stantons, so that the equation for the evolution 
of N = Ni + Nj is: 

I? =-4p 4 K)(^K) 2 -A7v 2 ). (7) 



If we assumed p(n c ) = cte, the solution of Eq. ([?]) 
would give N(n c ) ~ , the expression used 

in [HI , as a first order approach. But our cooling 
procedure modifies instanton's size (See figure ||), 
in a way that we phcnomenologically parametrize 
as: 



P {n c ) = p(0)(l + <zln(l + n c )) . (8) 

We will include (||) in equation (J?]), with p(0) the 
extrapolated radius at the thermalised situation 
and a a constant to determine. 

After performing a combined fit of our lattice 
results to the expressions (||) and the one coming 
from the integration of (Q), we can fix the initial 
values of the density, n(0) and the radius p(0), 
and the two constants that govern the evolution, 
A and a. An indirect measure of the (j4 2 ), . 
condensate is obtained. 

Nevertheless, the result of the extrapolation is 
highly dependent on the value of p(0), which due 
to its logarithmic behaviour is hardly reliable. We 
therefore prefer the value of {A 2 Ins ) at the maxi- 
mum of figure [| as a crude estimation of the in- 
stantonic contribution to (^4 2 )- 

Having a slight dependence on the lattice spac- 
ing, the result given here for this quantity has to 
be understood as a first estimation, that would 
require for a further analysis ]T^ |. 

4. Conclussion. 

The result of the combined fit gives a value 
of the instantonic contribution to {A? Ins ) ~ 
0.4GeV 2 , however, as reasoned before, we will 
give the value at the maximum as a more reli- 
able, although crude, estimation of the instan- 
tonic contribution to the {A 2 ) condensate. So 
our final result will be (Aj ns ) = 1.10(10) GeV 2 
for the simulation at lattice spacing 0.101 fm, 
and 1.37(16) GeV 2 for the one at 0.074 fm. 

This semiclassical evaluation of {A 2 ), which 
does not run with the scale, is difficult to re- 
late to that appearing in the O.P.E. expansion, 
which does depend on the renormalisation scheme 
and scale. The typical scale of instantons is 
p^ 1 ~ 0.7GeV. Unluckly it is not possible to run 
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Figure 3. Results of the combined fit for the instan- 
ton density and radius as a function of the number of 
cooling sweeps for a 24 4 lattice at f3=6.0. 



the (Aq p e ) to such a low energy, where per- 
tubative QCD is not valid. The lowest reacheble 
energy scale is 2.6GeV pjl^l, where the (A 2 ') con- 
densate takes the value; 



(^, P . s .(2.6GeV)) 



1.4(3)(3)GeV 2 



(9) 



the hrst error coming from the OPE determina- 
tion of the condensate renormalised at 10 GeV, 
and the second from higher orders in the running. 

Keeping in mind the level of uncertainty of 
these calculations, we can nevertheless claim a 
rather encouraging agreement between the in- 



stantonic contribution to the condensate and the 
one computed from the running of the Green 
Functions. 
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